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ABSTRACT 

Young A and B stars, the so-called Herbig Ae/Be stars (HAeBe), are surrounded by an active accretion disk and drive outflows. We 
study the jet HH 409, which is launched from the HAeBe star HD 163296, using new and archival observations from Chandra and 
HST/S7IS. In X-rays we can show that the central source is not significantly extended. The approaching jet, but not the counter-jet, 
is detected in Lya. In addition, there is red-shifted hya emission extended in the same direction as the jet, that is also absent in 
the counter-jet. We can rule out an accretion or disk-wind origin for this feature. In the optical we find the knots B and B2 in the 
counter-jet. Knot B has been observed previously, so we can derive its proper motion of (0.37 ± 0.01)" yr" 1 . Its electron density is 
3000 cirT 3 , thus the cooling time scale is a few months only, so the knot needs to be reheated continuously. The shock speed derived 
from models of Ha and forbidden emission lines (FELs) decreased from 50 km s -1 in 1999 to 30 km s -1 in 201 1 because the shock 
front loses energy as it travels along the jet. Knot B2 is observed at a similar position in 2011 as knot B was in 1999, but shows a 
lower ionization fraction and higher mass loss rate, proving variations in the jet launching conditions. 
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£^1. Introduction 

^ Stars form when molecular clouds fragment and contract to 
+jf proto-stars. Mass accretion onto those stellar cores proceeds 

jrt via an accretion disk, while the surrounding envelope disperses, 
i— i Eventually, the stars become optically visible; at this stage the 
low-mass population is called classical T Tauri stars (CTTS), the 
T 7~^ A and B star progenitors are Herbig Ae/Be stars (HAeBe). The 
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angular momentum accreted from the disk can be removed by 



outflows, which come in different types: There are slow, wide- 
angle winds, often seen in molecular lines such as H2, faster 
fV~) winds, showing up in optical forbidden emission lines (FELs) 
• such as [Oi], and highly collimated jets, which often reach ve- 
£j locities up to 400 km s 1 (Eisloffel & Mundt 1998). In very 
young objects these jets can be traced out to a few parsecs and are 
still well collimated. Internal working surfaces and bow shocks 
• • are visible in emission lines as knots, the so-called Herbig-Haro 
^ (HH) objects. The evidence for X-rays from outflows is rela- 
tively recent, starting with HH 2 (Pravdo et al. 2001; Schneider 
_j et al. 2012) and HH 154 (Bally et al. 2003; Favata et al. 2006; 
& Schneider et al. 2011). In both cases the jet driving sources are 
deeply embedded. 

To date, we know four CTTS or HAeBe with spatially ex- 
tended X-ray emission: The CTTS DG Tau (Giidel et al. 2005, 
2008, 201 1; Schneider & Schmitt 2008; Gunther et al. 2009) and 
RY Tau (Skinner et al. 201 1), Z CMa, a young star in an FUOr 
outburst phase (Stelzer et al. 2009) and the HAeBe HD 163296 
(Swartz et al. 2005). 

Jets from CTTS consist of layers with different velocities, 
where the innermost and most collimated component is the 



X 



* Figures 1-2 are only available in electronic form via 
http://www.edpsciences.org 



fastest one (Bacciotti et al. 2000) with velocities of a few hun- 
dred km s _1 . The outer layers are slower and less collimated. 
Often there is an asymmetry between the jet and the counter jet. 

Little is known about the physical mechanism, which drives 
the outflows. Different theoretical models of stellar winds (Kwan 
& Tademaru 1988; Matt & Pudritz 2005), X-winds (Shu et al. 
1994) and disk winds (Blandford & Payne 1982; Anderson et al. 
2005) have been proposed. Ultimately, the jet launching must be 
powered from the gravitational energy released in the accretion 
process. This is supported by the observation that the outflow 
rate is roughly one tenth of the accretion rate (Cabrit et al. 1990; 
Hartigan et al. 1995; Coffey et al. 2008), but it is unclear how the 
energy is converted. 

All models of jet launching and collimation rely on mag- 
netic fields of some kind. However, jets are observed not only 
from solar-mass CTTS, but also from HAeBes. In contrast to 
CTTS with magnetic fields typically in the kG range, HAeBes 
are not expected to posses outer convection zones; thus, they 
cannot drive a solar-like dynamo, although primordial magnetic 
fields might be present. Often no or only weak fields can be ob- 
served (Hubrig et al. 2007; Wade et al. 2007); the formal limit 
for HD 163296 is -25 ± 27 G (Hubrig et al. 2007). 

In this article we present new results on the time-evolution 
of the jet from the HAeBe HD 163296. 

In Sect. 2 we summarize previous work done on HD 163296 
and its jet HH 409. In Sect. 3 we give details on observations 
and data reduction. The results from the individual observations 
are given in Sect. 4, followed by a discussion (Sect. 5). We end 
with a summary in Sect. 6. 
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2. HD 163296 and its jet HH 409 

HD 163296 is surrounded by an accretion disk with a mass ac- 
cretion rate of M acc « 10~ 7 M G yr _1 . This disk shows a position 
angle of (140 ± 5)° and an inclinaton of i = (60 ± 5)° (Garcia 
Lopez et al. 2006) 1 . Perpendicular to it lies a chain of nebulosi- 
ties (position angle (42.5 ± 3.5)° Grady et al. 2000), called HH 
409. This jet is detected both by coronagraphic imaging (Grady 
et al. 2000; Wassell et al. 2006) and in HST/STIS long slit spec- 
tra (Devine et al. 2000). 

The distance to HD 163296 is 119^[q pc (van Leeuwen 
2007). The presence of a late-type companion can be ruled-out 
with a high degree of confidence from the HST imaging, the UV 
and optical spectra (Swartz et al. 2005, and references therein) 
and also more recent high-resolution optical spectroscopy is 
fully consistent with HD 163296 being single (Montesinos et al. 
2009). 

The jet can be traced down to 0'.'06 (7.1 AU) in the HST/STIS 
image and out to « 4000 AU in narrow band optical imaging 
(Wassell et al. 2006). Jet and counter jet are asymmetric. They 
have a different number of knots and different opening angles, 
but surprisingly the mass loss rate agrees better than a factor of 2 
on both sides (M jet « 10~ 8 M Q yr~' « 0.1 x M acc ). The opening 
angle is about 2° for the jet; the initial value for the counter jet 
is 5°, but it is more collimated further out (Wassell et al. 2006). 
The total space velocities, calculated from proper motion and 
radial velocity, are 360 km s _1 and 260 km s _1 , respectively. 
The knots of the jets radiate in Lya, Ha and optical FELs, e.g. 
[O i], [S n] and [N n]. Line ratios of the last three ions are density 
and temperature sensitive and thus give estimates for the electron 
densities in the HH objects of the order 10 3 cm 3 . Wassell et al. 
(2006) show structural differences between the knots: In knot 
C the Ha and [S n] emission are coincident, but in knot A Ha 
precedes the [S n] emission. We use the naming convention of 
Wassell et al. (2006) throughout this work to refer to individual 
knots. 

HH 409 seems similar to jets from CTTS in mass loss rate, 
electron excitation and opening angle. However, these quanti- 
ties are all measured from a jet component which moves at an 
intermediate velocity of 200-300 km s _1 . There are two indica- 
tions that the fastest component of HH 409 is more energetic. 
Both come from X-ray observations: First, Swartz et al. (2005) 
detected X-ray emission (luminosity 3 x 10 27 erg s _1 ) close to 
knot A at 7 "(=830 AU) from the star. Second, the X-ray spec- 
trum of the star itself shows the UV-sensitive O vn triplet, whose 
line ratios are not compatible with emission close to the stel- 
lar surface. This suggests X-ray emission from jet collimation 
shocks (within a few AU of the star). Using the observed X-ray 
temperature, Gunther & Schmitt (2009) estimate a flow veloc- 
ity of 500 km s _1 for the fastest jet component. Yet, the optical 
FELs and the Lya/Ha ratio indicate shock speeds of the order 
80-90 km s (Wassell et al. 2006), which is far too slow to pro- 
duce any X-ray emission. 

3. Observations and data reduction 

Table 1 gives a summary of the observations we used. 
3.1. Chandra 

The X-ray observations in table 1 were performed with Chan- 
dra/AClS-S in the VFAINT mode. HD 163296 is placed on the 

1 All position angles in this work are measured North to East. 



back-illuminated chip S3, because this chip has the best sensi- 
tivity for soft X-rays. This ACIS chip has a spatial resolution 
of 0.492" per pixel and an intrinsic energy resolution of 110- 
140 eV for energies below 2 keV. We reprocessed all Chandra 
data with CIAO 4.4 (Fruscione et al. 2006). Light curves and 
spectra for the central source are extracted using standard CIAO 
scripts. The background was determined from a large, source- 
free region on the same chip. Spectral fitting was done with the 
SHERPA fitting tool (Doe et al. 2007). 

The ACIS point-spread function (PSF) is under-sampled by 
the pixel size. To avoid an aliasing on the pixel grid two tech- 
niques can be used. Either the photon positions are randomized 
with offsets up to half a pixel or their positions are corrected with 
an energy-dependent sub-pixel event repositioning (EDSER, Li 
et al. 2004). The first method leads to a well-characterized, but 
artificially broadened PSF; the second method is less well char- 
acterized and was implemented in CIAO only recently, but it 
increases the resolution of the resulting images directly. All 
images shown here are processed with the EDSER algorithm, 
but we verified that images with randomized pixel positions give 
compatible results. 

The background light curves for both observations are flat 
confirming that no flaring of the ambient spacecraft proton im- 
pact rate was present. 

3.2. HST/STIS 

The jet of HD 163296, HH 409, has been observed with HST in 
different instrument configurations. Here, we analyze archival 
and new STIS long-slit spectroscopy, where the slit is aligned 
with the jet axis. The observations cover two different wave- 
length ranges: The G14QM grating is used for FUV spectroscopy, 
most importantly the hydrogen Lya line, and the G75QL and 
G750M gratings provide information on optical FELs with differ- 
ent wavelength resolutions. Table 1 gives details of the datasets 
used. In most cases the slit is narrow (52"x0'.'2), thus the exact 
position angle is important, only observation O66Q02010 was 
taken with a wider slit (0'.'5 versus 0"2). Grady et al. (2000) 
determine a jet orientation of 42.5 ± 3.5°. 

HD 163296 is optically bright (m v = 6.9 mag). To avoid 
saturation of the CCD the central star is placed behind the coro- 
nagraphic wedge (width 0.5", aperture 52X® . 2F1) in all optical 
spectra. The same setup was chosen for our new FUV spectrum 
as well. 

3.2.1. New optical observations 

Pipeline processed STIS data were downloaded from MAST. A 
full detector image is shown in Fig. 1. Optical CCD images 
are corrected for overscan and bias, then dark-subtracted, flat- 
fielded, rectified and wavelength calibrated to a heliocentric ref- 
erence frame. The spectra span the wavelength range from ap- 
proximately 6300-6865 A and the resolving power is R « 5000. 

Despite the coronagraphic wedge, there is a significant stel- 
lar contribution in the optical data (Airy rings). Below we show 
position-velocity diagrams (PVDs) of the strongest lines in our 
data. To remove the stellar contribution we define two wave- 
length intervals to the left and to the right of the emission line of 
interest (The specific wavelength regions chosen are given with 
the PVD). We then estimate the stellar contribution in the re- 
gion of the PVD by a linear interpolation between the left and 
the right wavelength region. This procedure works because the 
wavelength intervals of interest are very narrow. However, low 
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Table 1. Log of observations 



Facility 


Instrument 


Id 


exposure 


date 






P.A. 


central wavelength 


aperture 


Chandra 


ACIS-S 


3733 


20 ks 


2003- 


-08- 


■10 








Chandra 


ACIS-S 


12359 


46 ks 


2011 


-02- 


■09 








HST 


STIS/G750L 


O5FO22010 


1140 s 


1999- 


-08- 


■09 


43.06° 


7751 A 


52X0.2F1 


HST 


STIS/G140M 


O57Z03010 


2230 s 


1999- 


-09- 


■30 


45.15° 


1218 A 


52X0.2 


HST 


STIS/G140M 


O66Q02010 


2282 s 


2000- 


-07- 


■22 


46.25° 


1218 A 


52X0.5 


HST 


STIS/G140M 


O66Q03010 


7707 s 


2000- 


-07- 


■21 


46.50° 


1218 A 


52X0.2 


HST 


STIS/G750M 


OBI801010 


2275 s 


2011- 


-07- 


■21 


44.06° 


6581 A 


52X0.2F1 


HST 


STIS/G140M 


OBI801020 


2700 s 


2011- 


-07- 


■21 


44.00° 


1218 A 


52X0.2F1 



surface brightness features in the jet cannot be identified close 
to the central object due to the low contrast to the bright Airy 
rings. Furthermore, the strong, presumably stellar, Ha profile 
scatters into the outer jet regions and appears as extra emission 
in the continuum-subtracted position-velocity diagrams. A com- 
parison with spatial profiles close to the Ha wavelength reveals 
that all excess flux is compatible with scattering from the central 
position. 

3.2.2. New FUV observations 

Due to a change in the Science Data Formatter (SDF) on HST 
not all events in the new G140M exposure were transferred from 
the SDF to the on-board solid-state recorder (SSR) for downlink. 
This leads to a gap of 138.5 s in the exposure. Once the damaged 
period is removed from the telemetry stream, the resulting raw 
data is processed with the standard calstis pipeline in the same 
way as the raw data from all other observations. A full detector 
image is shown in Fig. 2. 

CalSTIS divides the image by a flat-field, then it is rectified 
and wavelength calibrated. The data covers approximately the 
range 1190-1245 A. We use a binning of 4 high-resolution pix- 
els, which provides a wavelength resolution of 0.2 A for a point 
source. For comparison the effective resolution of emission that 
fills the slit uniformly, would be w 80 km s _1 . 

While the overall dark current is low, the top left hand corner 
of the FUV-MAMA detector shows an enhanced glow that is 
intrinsic to the micro-channel plate array. The strength of the 
glow is time variable, but it can be removed by scaling archival 
dark exposures. However, these have a limited signal-to-noise 
ratio and thus the subtraction would increase the noise in the 
residual image. Since the glow changes on much larger spatial 
scales than the features of the jet, we fit the local background 
around the knots of the jet instead of subtracting it globally. 

The FUV images are dominated by strong Lya airglow emis- 
sion. For each position we determine the centroid of the airglow 
line. We build an average line profile from regions close to the 
detector edges and far from the known jet emission. At each 
position we subtract this profile from the image. The accuracy 
of this procedure is limited by the statistical noise in the air- 
glow line. We estimate that only features with a surface lumi- 
nosity of > 10% of the observed airglow line can be detected. 
All residuals are smaller than this value. The data is taken in 
TIME-TAG mode, and as a cross-check we limited the analysis to 
times where the count rate in the airglow is less than a third of 
the maximum rate. We conclude that no astrophysical signal is 
found at the position of the airglow line. 

In the FUV the Airy rings are much less severe. The gener- 
ally lower count rates mean that the interpolation is not as well 
defined and the procedure described above would introduce sig- 



nificant noise to the PVDs. Therefore, no subtraction of the con- 
tribution from the unresolved source was attempted for the FUV 
exposures. 

3.2.3. Archival data 

Pipeline processed STIS data were downloaded from MAST and 
treated in the same way as the new FUV or optical data. While 
all G140M observations are useful to analyze material close to 
the central source, we use only the observations from 1999 and 
201 1 when we look at the knots of the jet, which are located a 
few arcsec from the star because Wassell et al. (2006) show that 
the position angle of 46.5° in the deeper observation O66Q03010 
is too large and the slit misses the apex of the jet shocks. Obser- 
vation O66Q02010 is the only dataset taken with a wider slit. 
While knots are visible, their flux is uncertain because they par- 
tially overlap with the strong Lya airglow. 

3.2.4. Data reduction 

We further processed the data with custom python routines. Hot 
or negative pixels are marked by hand and replaced with the av- 
erage value of the surrounding pixels. The jet emission is weak 
both in the FUV and in the optical. Consequently, the signal- 
to-noise ratio (SNR) for a single pixel is typically 2-4 in the jet. 
However, all features analyzed stretch over a larger detector area. 
We present the data binned in the spatial direction over 4 pixels 
for the FUV and over 6 pixels for [O i] 6300A and the [S n] lines 
at 6716A and 6731 A. When we fit fluxes and positions for knots 
in the jet, we achieve even better SNR because they stretch over 
a range of velocities as well. Uncertainties given in the tables 
are lcr confidence intervals throughout the article. 

4. Results 

4. 1 . X-ray emission 

The first Chandra observation of HD 163296 was already pre- 
sented by Swartz et al. (2005) but we analyze it here again for 
comparison. HD 163296 was also observed with XMM-Newton 
by Gunther & Schmitt (2009). 

4.1.1. The central star 

The light curves and hardness ratios for the archival observations 
are moderately variable and the new observation is fully compat- 
ible with those findings. The XMM-Newton spectra have a higher 
signal-to-noise ratio than the Chandra data because the effective 
area of XMM-Newton is larger, the exposure time is longer and 
XMM-Newton/RGS delivers a high-resolution grating spectrum. 
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Table 2. Best-fit model parameters (l<x confidence interval). 



component 


soft 


medium 


hard 


Parameters fixed for all components' 1 


N H 


[10 2l) cm" 2 ] 


7 


7 


7 


P 


b 


3.7 




J. 1 





b 


0.7 


0.7 


0.7 


Ne 


b 


1.2 


1.2 


1.2 


Mg 


b 


2.3 


2.3 


2.3 


Si 


b 


2.8 


2.8 


2.8 


Fe 


b 


1.6 


1.6 


1.6 




Chandra 2003 (red. x L 


= 0.91) 




kr 

VEM 


[keV] 

[10 52 cm" 3 ] 


19+ - 23 

1 3 +0 - 5 


60 +0(b 

"•""-0.04 

9 +0 - 2 

"•^-0.2 


2 +1 ° 

4 +M 
"•^-o.i 


XMM-Newton 2007 c 


kr 

VEM 


[keV] 

[10 52 cm" 3 ] 


"•^-0.01 
o 9+0.5 
^ -0.3 


51+ 01 

u - Ji -0.02 

1 l +0 ' 3 
1,1 -OA 


2 7 +10 
5^ : ' 

u -0.2 




Chandra 201 1 (red. \ l 


= 0.85) 




kr 


[keV] 

[10 52 cm" 3 ] 


l6 +om 

Q oi 

1 3 +0 - 2 


60 +0 - 2 

"•""-0.02 

1 o +0 - 6 


2 4 +() - 5 

4 +M 
"• -0.1 



(a) from the XMM-Newton fit. See Gunther & Schmitt (2009) 

for details. (b) abundances relative to the solar values of 
Grevesse & Sauval (1998) (c) All values taken from Gunther & 
Schmitt (2009). 



We use the same spectral model for HD 163296 as in (Gunther 
& Schmitt 2009) (one phabs absorber and three vapec, opti- 
cally thin, collisionally excited, thermal emission components 
with non-solar abundances) and keep the abundances fixed at 
the values found in the XMM-Newton observation. Elements like 
Mg, Fe, Si and C with a low first-ionization potential (FIP) are 
enhanced compared to O and Ne which have a higher FIP. This 
is a common effect in inactive stars. Also, we keep the absorb- 
ing column density fixed at 7Y H = 7 x 10 20 cirr 2 , which is in full 
agreement with the optical reddening (Gunther & Schmitt 2009). 

Table 2 lists the temperature and volume emission measures 
for the three temperature components. Both Chandra observa- 
tions are refit here, the reduced x 2 values of both fits are below 
1. Statistically acceptable fits are possible with two temperature 
components only, but we keep the model of Gunther & Schmitt 
(2009) to facilitate a comparison. The fit values for the XMM- 
Newton observation are taken directly from that work (after cor- 
rection for the revised distance we use in this article). The tem- 
peratures of all three components are very similar and consistent 
within the errors over all three observations. The same is true for 
the volume emission measure (VEM) with the exception of the 
soft component, which seems marginally stronger in the XMM- 
Newton observation. 

The light curves show little variability during the obser- 
vations and the comparison of the spectra and their emission 
measures over nearly one decade shows that the X-ray emis- 
sion of the HAeBe HD 163296 is stable. The observed (ab- 
sorbed) energy flux in the 2011 observations is (2.4 ± 0.2) X 
10" 13 erg s _1 cirr 2 (0.3-5.0 keV). From this we derive an in- 
trinsic \ogL x = 29.8 (in ergs s -1 ), which is almost identical to 
the values of previous observations (Gunther & Schmitt 2009). 

4.1.2. The base of the jet 

Based on a line ratio in the Ovu He-like triplet that is sensitive 
to the ambient UV field Gunther & Schmitt (2009) conclude that 



the soft X-ray emission must be emitted at least 0.8 R* above the 
stellar surface. 

In the CTTS DG Tau Schneider & Schmitt (2008) directly 
measure an offset of about 40 AU between the centroids of the 
soft and the hard component of the central source. The hard com- 
ponent is presumably stellar emission, while the soft component 
might originate in a jet shock close to the star. This detection is 
possible in DG Tau, because the high column density absorbs all 
soft X-rays from the star itself. We search for an offset using the 
same technique in HD 163296 but find no hint of extension. A 
larger physical separation would be required in HD 163296 than 
in DG Tau, because the absorbing column density is very low 
and so any stellar plasma would also contribute to the soft band 
flux. 

Additionally, we simulate the Chandra point-spread function 
with SAOTrace Version 2.0.1 (Carter et al. 2003). This tool per- 
forms Monte-Carlo ray-tracing through the Chandra optics. In 
the current version it includes the telescope dither pattern. From 
the ray-trace we calculate the detector response with MARX 
(Wise et al. 1997), version 5.0. The combination of SAOTrace 
and MARX can properly simulate the effect of the EDSER sub- 
pixel event repositioning. We find an extension beyond the sim- 
ulated PSF in the range 0.5-0.8". However, a deconvolution of 
the observed image with the simulated PSF using the CIAO tool 
arestore shows that the extension in both observations is per- 
pendicular to the jet direction. It matches the PSF asymmetry re- 
cently detected in Chandra/HRC observations (Juda & Karovska 
2010) in both the distance to the center of the PSF and in the rel- 
ative position in detector coordinates. Thus, the central source it- 
self is apparently not extended. In comparison to DG Tau, where 
the jet is detected within 30 AU of the central star using image 
deconvolution (Giidel et al., in prep.) we conclude that almost 
all soft X-ray emission from HD 163296 is contained within the 
inner 30 AU. 

4.1 .3. The X-ray knots of the jet 

Figure 3 shows the region around HD 163296 in X-rays. The po- 
sition angle of the jet is well known from optical images. Thus, 
we can select the jet region a-priori and do not need to rely on a 
source detection in the Chandra images. We select events with 
energies in the range 0.2-2.0 keV. Fig. 4 shows a projection of all 
counts within 1" of the jet axis. The black lines show the count 
distribution along the jet. For comparison, the same method was 
employed for a set of angles rotated by 20°, 40°, 60°, 80°, 100°, 
120°, 140° and 160° with respect to the jet (red lines). 

In the data from 2003 the black line is above the compar- 
ison lines in the range -3 to -4" (negative numbers are for the 
approaching south-west jet). To estimate the significance, we 
determine an average number of 1.75 counts between -4.5" and 
-2.5" from the wings of the PSF, while 5 counts are observed in 
this range along the jet direction. In a Poisson distribution the 
chance to observe 5 counts when the expectation value is 1 .75 is 
3%. Thus, this feature is 97% significant. Note, that this does 
not contradict the results of Sect. 4.1.2, because this feature is 
two orders of magnitude less luminous than the central star and 
thus it does not show up in an image deconvolution. This emis- 
sion corresponds to the extension of the source seen by Swartz 
et al. (2005). A similar feature in 201 1 is seen in the black line 
but also in directions other than the jet (red lines in the figure), 
thus it is not clear if this feature is present in 201 1. 

In the projection the strongest features beyond the stellar PSF 
are two peaks at -7.2" and -10.6" in 2003 and 201 1, respectively. 
Swartz et al. (2005) derive a probability of only 0.04 that the knot 
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2003 - ObsID 3733 




22.0s 21.5s 21.0s 20.5s 17h56m20.0s 
RA (J2000) 



2011 - ObsID 12359 



Q 




22.0s 21.5s 21.0s 20.5s 17h56m20.0s 
RA (J2000) 



Fig. 3. Chondral ACIS image of HD 163296 
and its jet HH 409 in the energy range 0.2- 
5.0 keV binned up by a factor of two (Each 
block is 2 pixel « 1" on the side.). The black 
box shows the position of the optical jet. The 
half-width of the box is comparable to the PSF 
(1"). The long axis of the box passes through 
the centroid of the X-ray emission and has the 
position angle 42°. Circles indicate locations 
of possible emission knots. The radii are 2.5". 
See the electronic edition of the Journal for a 
color version of this figure. 



in 2003 (black circle in the left panel of Fig. 3) is a chance fluctu- 
ation by calculating the number of circles of this size at any posi- 
tion on the detector that contain 5 or more photons. We confirm 
this number using Poisson statistics of the measured background 
in a near-by source free region (0.008 and 0.012 cts pixel -1 in 
2003 and 2011, respectively, in the energy range 0.2-2.0 keV). 
Thus, we expect 0.65 and 0.97 background events in a circle with 
radius 2.5". According to Poisson statistics the chance to find 5 
or more photons in any one circle in 2003 is therefore < 0.1% 
and the chance is 7% to find 3 or more photons in 2011. To es- 
timate the significance of the detection, we need to consider the 
number of circles searched: We analyze a region up to 25" from 
the central star in both jet directions, which can be covered by ten 
circles with radius 2.5". (We would recognize features that fall 
between two of those circles, thus more than ten positions are 
searched. However, only the number of independent positions 
is important.) In 2003, the probability to find no circle with 5 
counts or more is 99%, which we can take as the significance of 
the detection of an X-ray feature in the jet; the significance of 
an independent detection of the X-ray feature in the jet in 2011 
is only 44%. A potential knot with three photons in 201 1 is not 
significant by itself. However, if we assume the X-ray feature 
in 2003 is real and then propagate it with the proper motion of 
the optical knot A, which shows 0"49 + 0"02 yr -1 (Wassell et al. 
2006), then this defines the search area exactly (circle in Fig. 3, 
right panel) and the chance to find three of more photons in a 
single circle of this size is < 7%. Thus, we can say that as- 
suming the knot in 2003 is real and has the same proper motion 
as the optical knots, we detect an X-ray knot in 201 1 with 93% 
significance. 

We note, that the half-width of the detection region used by 
Swartz et al. (2005) is 2.5", much larger than expected from the 
width of the optical knots and the Chandra PSF. Fig 4 shows the 
projection of the counts on the jet axis with an extraction region 
matched to the PSF. Several features in the comparison sample 
(projected on lines which do not coincide with the jet axis) have 
similar fluxes as the tentative knots along the jet. However, in 
HH 2 Schneider et al. (2012) find that the X-ray emitting region 
is larger than the optical knots. If the same is true here, then it is 
not surprising, that the knots are only detected if we use a large 
extraction region. 

The effective area of the Chandra/ ACIS detector degraded 
between 2003 and 201 1 (about 5% at 1 keV, 20% at 0.8 keV and 
40% at 0.5 keV), but the observation taken in 201 1 is more than 
twice as long. Taking the three photons at -10.6" as an upper 
limit to the flux of the tentative X-ray knot, the apparent lumi- 
nosity decreased by a factor three compared to 2003. Although 
the count number is low the hypothesis that the photon flux is 
the same in both observations can be rejected on the 95% level. 
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Fig. 4. Black line: Projection of Chandra/ ACIS counts onto the jet 
axes in a rectangle with half-width 1". red/gray lines: projection onto 
a set of lines through HD 163296, but with different angles. All fluxes 
are smoothed. Distances are positive for the NE jet and negative for the 
SW jet. See the electronic edition of the Journal for a color version of 
this figure. 



In summary, we can say that in 2003 an X-ray knot is de- 
tected with 95% significance, confirming Swartz et al. (2005). 
Under the assumption that this knot moves with a proper mo- 
tion similar to the optical knots (as observed in HH 2 Schneider 
et al. 2012), this feature is detected again in 201 1, but with lower 
luminosity. 



4.1 .4. Summary of X-ray results 

HD 163296 itself is very stable over the last decade. Lightcurves 
and hardness ratios in all observations are nearly flat and the tem- 
peratures and emission measures of all spectral fits agree. From 
the O vii line ratio Gunther & Schmitt (2009) conclude, that the 
soft X-ray emission must originate above the stellar surface, but 
we do not see an offset between the soft and the hard centroid, 
setting an upper limit on the separation of 30 AU. However, very 
weak emission features (knots of the jet) are seen on the edge 
of the PSF at -4" along the jet. Furthermore, there is a weak 
feature at -7.2" in 2003 and a tentative feature at -10.6" in 201 1 
(table 3). Both of these features do not have an optical or UV 
counterpart (see below). 
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Table 3. Position of knots from Gaussian fit 



knot 


tracer 


epoch 


position 


X-ray knot 


V ran 


70fn £>Q 


— / .Z ± 1 ,u 


X-ray knot 


X-ray 


2011.11 


-10.6±0.3" a 




Lya 


2011.64 


* -2" 




Lya 


2011.64 


-4.0 ± 0.2" 


knot B 


[Sn] 6731 A 


1999.69 


4.3 ±0.1" 


knot B 


[Sn] 6731 A 


2011.64 


8.7 ±0.1" 


knot B2 


[Sn] 6731 A 


2011.64 


3.3 ±0.1" 



(a) The uncertainty is calculated from the variance of the photon 
position for all photons attributed to the knot. 

4.2. UV emission 

HH 409 was observed at UV wavelengths with STIS previously 
in 1999 and in 2000. The most prominent line in the FUV is Lya. 
Its central part is contaminated by ISM absorption and airglow, 
which fills the aperture. Thus, only data more than 150 km s _1 
from the line center can be used. Devine et al. (2000) found blue- 
shifted Lya emission in the jet with a velocity around 380 km s _1 
at 0.06" (7.1 AU), that decreased to 335 km s" 1 at 6" (714 AU). 
Additionally, they could identify three knots (Knot A in the jet 
at -7.7" and knot B and C in the counter-jet at 4.5" and 8.6", 
respectively). 

Figure 5 (left) shows the position-velocity diagram (PVD) of 
Lya for the new observation. Despite the slightly longer expo- 
sure time the data quality of the new observation does not reach 
the level of the observation in 1999. 

In the jet itself there are two knot features at about -2" and 
-4". At the given SNR ratio they are well described by a Gaus- 
sian. Fitted values for the position are given in table 3, line fluxes 
in table 4. Figure 6 shows spectra extracted at the position of the 
two knots. This is the same range of velocities and speeds ob- 
served by Devine et al. (2000) and indicates that the flow speed 
of the jet has not changed considerably in the last decade. 

The Si in line at 1206.5 A is too weak to be fitted indepen- 
dently. Thus, we fix the line shift and width at the values found 
for Lya. For the outer knot we can only give an upper limit, for 
the inner one we formally find a detection (table 4). However, 
we caution that the flux is weak and the position of the inner 
knot coincides with a region of enhanced detector background 
in Si m. Thus, the flux might be overestimated. 

In principle the models of Hartigan et al. (1987) constrain 
the shock velocity given the Lya/Si in flux ratio. These model are 
based on simulations of radiative shocks. The pre-shock material 
is assumed to be in ionization equilibrium with the post-shock 
material. This describes a scenario where the Lyman radiation 
from the shock front is the main ionization source for the pre- 
shock gas. At a distance of a few hundred AU from the central 
star this is a reasonable assumption. 

In all these models the flux ratio Lya/Si m> 10. Thus, they 
are not constrained by the upper limit in the outer knot and they 
cannot explain the tentative Si m detection in the inner knot. It is 
possible that scatter from the central star contributes to the inner 
emission similarly to the spatially extended Lya emission, which 
is seen red-shifted on the side of the approaching jet and which 
possibly also contributes to the blue-shifted emission. 

The observed surface brightness of the Lya emission is re- 
duced compared to 1999. This can be seen in table 4 that shows 
the fluxes integrated over the visible knot surface. Either the jet 
is dimmer in 201 1 or it has broadened considerably, so that only 
part of the total luminosity is contained in the aperture. The later 
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Fig. 7. Spatial profile on the red (250 to 450 km s~') side of the 
Lya line, the Sim line and of a continuum region. All line profiles 
and continuum profiles are normalized (see text for details). Shaded 
regions indicate the statistical uncertainty for the lines. The uncertainty 
on the continuum profile is negligible because it is extracted from a 
much larger region. 



explanation seems unlikely, because the jet was centrally bright 
in previous observations. Also, in a quasi-static scenario where 
gas travels along the jet but the overall shape of the jet stays 
constant with time this is not expected. 

The line fluxes in table 3 are not dereddened. Devine et al. 
(2000) found E(B - V) = 0.015, which leads to flux correction 
of 15% (Wassell et al. 2006). 

Given the proper motion of the old knots A and C (Wassell 
et al. 2006) they should be located at -13.6" and 12.7" during the 
new observation. No significant emission is observed at these 
locations, nor at the tentative X-ray knot at -10.6". 

After subtracting the spectrum of the unresolved source 
Devine et al. (2000) find red-shifted Lya emission to the south- 
west of the central source in the direction of the blue-shifted jet. 
Figure 7 shows the spatial profile of the emission on the red- 
shifted side of the Lya line in comparison with the continuum 
emission, which originates from the star. The continuum and 
lines are normalized at 0.04", where the continuum has the max- 
imum of its spatial profile in our spatial binning. The red-shifted 
Lya emission shows extended emission on the side of the ap- 
proaching jet in all observations. The blue-shifted side is also 
extended (see Fig. 5), but since the jet is seen on this side we 
cannot distinguish if this emission is entirely due to the jet or if 
there is an additional component similar to that shown in Fig. 7. 
In contrast, the Si m line is consistent with the continuum, which 
is not spatially extended. 

The two observations taken in 2000 are less than a day apart 
and the spatial profiles have a very similar shape. Thus, it is un- 
likely that the difference in flux between them is due to temporal 
variability. Rather, it is caused by the aperture. One observation 
used a slit width of of 0.2", the other 0.5". Both spatial profiles 
are very similar, so the bulk of the emission in the direction per- 
pendicular to the slit and thus also perpendicular to the jet, must 
be contained within 0.2". 

We searched for a similar feature in the Sim line at 
1206.51 A, but the projections of line and background are fully 
compatible here. 
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Fig. 5. PVD for Lya-, Ha and three opti- 
cal FELs. For comparison, the contourlines for 
[S n] 6731 A are overlaid in all panels, left: The 
Lya- line is contaminated by airglow close to 
the rest velocity, right: The optical lines show 
scattered light close to the central star despite 
the occulting bar. These areas are shaded. [O i] 
6300A is located close to the edge of the chip, 
so no data exists in the bottom left of the PVD. 
Contour lines are drawn for [S n] 6731 A at the 
level 7 x lCT 16 erg s _I A -1 cm 2 arcsec~ 2 . The 
optical PVDs are background subtracted (see 
text). The following wavelength intervals are 
used to fit the background: [Oi] 6300A: 6352- 
6357A and 6376-6382A; Ha: 6504-6524A and 
6604-6624A; [Sn] 6716A: 6705-6708A and 
6721-6724A; [Sn] 6731 A: 6721-6723A and 
6735-6738A. See the electronic edition of the 
Journal for a color version of this figure. 
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Fig. 6. Spectra of the Lyo- line in the two brightest knots. The inner part of the line is contaminated by airglow, additionally the knot closer to 
the source (left panel) shows some emission on the red side of the line which could be due to scattered light from the central star (see text). We fit 
two Gaussians for knot and background (dashed line). Both knots appear blue-shifted by about 350 km s -1 . 



4.3. Optical emission 

PVDs of [Oi] 6300A, Ha, and the [Sn] lines at 6716A and 
6731A are shown in Fig. 5 (right). Stray light is strongly re- 
duced due to the occulting wedge used in the observations, but 
the inner few arcseconds of the long-slit data are contaminated 
by Airy rings caused by refraction of the bright central star. We 
fit and remove the Airy rings, but the residual image shows a 
higher noise level at these positions. 

On first sight, the Ha emission seems extended in the direc- 
tion of the blue-shifted jet - not only for the jet velocity, but 
indeed for a wide range from -300 to +300 km s _I . A projection 
similar to Fig. 7 shows, that this apparent extension is indistin- 
guishable from the spatial profile of the continuum. The only 
reason why it shows up prominently here is that the Ha line is 
far brighter than the stellar continuum. 

Two knots can be identified in the counter-jet; no feature is 
seen in the jet itself (table 3). The strongest knot (which we call 
B2) peaks at 3.3", very close to the outer edge of the disk as seen 
in scattered light (Grady et al. 2000). Even without contamina- 
tion by the Airy rings, the disk would block any emission from 



the counter-jet further in, thus we cannot say if the observed po- 
sition of B2 represents the peak of the knot or if we merely ob- 
serve the end of a more extended structure. The second knot (B) 
is located at 8.7". B2 and B show a similar spacing as B and C 
in the old observation (Wassell et al. 2006). 

Using the extraction regions 2.6-3.7" and 8.1-9.0" we ex- 
tract spectra for B2 and B, respectively (Fig. 8). In addition to the 
lines shown in the PVD we include [N n] 6586A. Gaussian fits to 
the lines are given in table 4. Ha is described by two Gaussians. 
The table reports the component with a higher velocity shift, as 
the inner component is most likely due to scattered stellar emis- 
sion. Due to the lower spectral resolution Ha and [N n], and the 
two [S n] lines at 6716 and 6731 A were blended in the previous 
optical STIS spectra. We remeasured the line fluxes from 1999, 
but our values are consistent with Wassell et al. (2006) within 
the errors. 

The line fluxes in the table are not dereddened. The wave- 
length range of the optical observations is small, so a very similar 
correction applies to all those lines. The same E(B — V) = 0.015 
as above leads to a flux correction of 3-4% in the optical. This is 
much less than the statistical errors on the line fluxes. 
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Table 4. Fluxes in knots (lcr uncertainties) 



Line 


A 


flux (1999) flux (2011) 


V 


cr 




A 


10~ 16 erg s _1 cirT 2 


km s _1 


km s _1 


jet: -3 to -1 " 


Lya 


1216 


49- 


353^ 


40 +z 


Sim 


1207 


16 a 


=353 a 


=40 a 


jet: -5.1 to -3.3" 


Lya 


1216 


32 ±2 


331™ 


771 


Sim 


1207 


< 4 a 


=331 a 


=77 a 


counter-jet: B2 


[Oi] 


6300 


6.3 + 0.6 


130+10 




Ha 


6563 


8.4 +5J 


144+3 

-13 


50 ±4 


[Nn] 


6586 


2 9 +0>: ^ 




33_ + 6 3 


[Sn] 


6716 


6 4+^ 

U - -0.5 


132 ±4 


31 ±3 


[Sn] 


6731 


9.1+0.6 


130 ±3 


35 ±3 



counter-jet: B 



TOi] 6300 2.3+J-i 2.0+"^ 139! 3 36+' 6 

Ha 6563 21 ± 2.0 b 2.0 + ° : ? 132 ± 4 20.02+^ 



+| 119 + 7 25.7^ 
[Sn] 6716 .__ +1 . 3 2.2 + °| 156 + 7 24 + 7 



[Nn] 6586 4.8+};^ 1.4™ 119 ±7 25.7_ 0J 



[Sn] 6731 ' 20 - 2 - 3 3.4 ±04 153+. 3 28^ 



=5 ""-3 



Notes. (a) The Si m lines is too weak to fit without constraints. Thus, we fix v and <x at the values found for Lya at the same position. <b) In the 
G750L grating used in 1999, the Ha line is blended with [Nn] 6548A. This contribution is subtracted assuming the quantum mechanical ratio 
[Nii]6584A/[Nii]6548A = 3.0. 



In knot B2 in the counter-jet the velocities of all lines 
are roughly compatible with a systematic knot velocity of 
130 km s -1 . The forbidden lines all have a similar width, only the 
Ha line is significantly wider (table 4). In knot B the [S 11] lines 
seem slightly more red-shifted than the other observed lines, 
however, the errors given in the table are statistical only. In gen- 
eral, the velocities in the optical are smaller and narrower than 
in the Lya line. Wassell et al. (2006) already found a velocity 
asymmetry between jet and counter-jet but at much lower level 
and with larger errors (223 ± 81 vs. 171+89 km s _1 ). Appar- 
ently, the Lya emission in the jet traces a faster component than 
the FELs. In our observation the FUV features in the jet are not 
detected in the optical, thus a comparison of FUV and optical 
lines in the same knot is not possible. Velocity asymmetries are 
not unusual. The ratio of jet velocity and counter-jet velocity is 
> 1 .5 for about half of all CTTS with optically detected outflows. 
Those differences seem to be intrinsic to the launching process 
and cannot be explained by interactions with the circumstellar 
matter close to the central source (Hirth et al. 1994). 

4.3.1 . Proper motion and velocity of knot B 

Knot B was seen at 4.3" in 1999. For reasonable values of 
its proper motion it should be located a few arcsec further out 
in the new STIS observation. Indeed, a knot is located at 
8.7" and we identify this as knot B with a proper motion of 
0.37 + 0.01" yr _1 which agrees well with the proper motion of 
knot C of (0.34 ± 0.02)" yr" 1 . Thus, it seems that knots in the jet 
move consistently with similar velocities. 

We take the velocity along the line of sight as the mean of 
all observed line shifts for knot B (table 4) and derive a total 
space motion of 250 km s _1 for knot B and an inclination of 
56°. The proper motion measures the motion of the shock front, 
while the line shift shows the motion of the gas. Thus, the above 
calculation requires that the shock front is co-moving with the 



material. As shown below, the relative shock speed is low and 
thus this approximation is valid. The measured values agree with 
previous measurements of the jet speed and inclination. 



4.3.2. Physical conditions in knot B 

While the observed surface luminosity in knot B2 is about 2- 
3 times higher than in knot B, the line ratios are very similar, 
which indicates a similar density and temperature, although knot 
B is nearly twice as far from the central source as knot B2 (inte- 
grated fluxes are given in table 4). A decade ago, when knot B 
was located close to the current position of knot B2, its surface 
luminosity in Ha and [S n] was significantly higher compared 
with knot B today. The surface luminosity of knot B in the [S n] 
and [0 1] lines reduced by a factor of 2, in [N 11] by a factor of 4, 
but in Ha by an order of magnitude between 1999 and 201 1. 



4.3.3. Electron densities 

The ratio of the [S 11] lines at 6716A and 6731 A is density sen- 
sitive with only a very modest temperature dependence. The 
differential reddening between those two lines is negligible, be- 
cause HH 409 is only weakly reddened and the wavelength is 
very similar. Unfortunately, these lines are unresolved in the 
early STIS observations, so that we do not have a compari- 
son value for knot B from an earlier epoch. However, in 2004 
the density of knot C was measured from ground-based data 
to n e = 1400 ± 400 ctrT 3 (Wassell et al. 2006). With the ob- 
served line values from table 4 we find n e — 2200 + |° ) ° ) ° cirT 3 and 
n e — 3000+ cirT 3 , for knots B2 and B, respectively, using the 
CHIANTI database (Dere et al. 1997; Landi et al. 2012). Be- 
low we calculate the total density n from n e and the ionization 
fraction. 
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Fig. 8. Spectra for Ha and several optical FELs at the position of the 
two brightest knots. The main component is fitted with a Gaussian. 



4.3.4. Shock velocities 

We estimate the shock velocities and the ionization fraction in 
the knots using the models of Hartigan et al. (1994). They con- 
tain updates for some optical lines over the simulations by Har- 
tigan et al. (1987), which we used to interpret the Lya7Sini ra- 
tio. Hartigan et al. (1994) present a series of plots for different 
line ratios that provide a diagnostic of the ionization fraction and 
the shock velocity. The models are calculated for three different 
pre-shock densities and for different magnetic fields. Given the 
estimate of the electron density above, we use the models for a 
pre-shock density of 10 3 cm 3 . The magnetic field in the shock 
zone is degenerate with the shock speed because a field paral- 
lel to the shock front inhibits strong compression. We assume a 
weak field (< 0.1 mG) here. 

Table 5 shows the ionization fractions < / > (averaged over 
the emitting area) and the shock velocities v s hock from the Harti- 
gan et al. (1994) models. The shock velocities derived from dif- 
ferent line ratios agree well. Knot B and B2 in the most recent 
observation have a very similar shock speed around 30 km s _1 , 
less than the shock speed of knot B in 1999. The differences be- 
tween the ionization fractions are much larger. Again, the knot 
B in 1999 was more highly ionized than any knot in 201 1, how- 
ever, the recent values span a range of a factor of five. This could 



Line ratio 



B(1999) B(2011) B2(2011) 



Ionization fraction 



[Oi]/Ha 
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0.02 
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[S n] (671 6+673 l)/Hor 
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0.02 


0.04 


[Nn]/[Oi] 


0.3 


0.1 


0.1 


shock 


velocity [km s~ 
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[Oi]/Ha 


50 


30 


40 


[S n] (67 16+6731)/Ha- 


40 


30 


30 


[Nn]/[Oi] 


70 


30 


30 


[Sii]6716/[Sii]6731 




45 


40 



be due to an underestimate of the Ha flux in the fitting process. 
In knot B2 the fitted Gaussian has a larger velocity than all other 
lines. Thus, a larger fraction of the flux is described by the Gaus- 
sian of a stationary component, which we interpret as scattered 
light from the central source. On the other hand, the relative 
errors on the line fluxes for weak lines are large and the differ- 
ences in < / > in table 5 are only marginally significant. In the 
following we adopt a shock velocity of 50 km s _1 for knot B in 
1999 and 35 km s _1 for knot B and B2 in 201 1. The shock mod- 
els predict a compression < C > of gas (averaged over the [S n] 
emitting area) of a factor 50 and 25 for these shock velocities. 

4.3.5. Mass loss rates 

There are several different methods to estimate the mass flux in 
a jet. In this section we follow Hartigan et al. (1994). First, the 
mass flux can be estimated from the total flux in a line, e.g. [O i]. 
We calculate the total luminosity and assume that the plasma is 
emitting close to the peak formation temperature of this line and 
that the total oxygen abundance is close to solar. This calcula- 
tion misses flux components that are located outside the aperture, 
ionized, too cool to emit [O i] 6300A or have densities above the 
critical density for [O i], e.g. in a dense and fast component very 
close to the jet axis. While we know little about temperature 
components outside of the formation range of FELs, the den- 
sity is most likely below the critical density (10 6 cm 4 ) given the 
estimates above. We use 
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which is eqn. 10 from Hartigan et al. (1994). Alternatively, a 
very similar calculation can be done using the [S n] lines. If they 
are in the high-density limit (> 10 4 cm -3 ), then the formula does 
not depend on the density in the knot (Hartigan et al. 1995, eqn 
A. 10): 
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4.5 x 10" 9 



sky 



10 16 cm 



-t-6731 

10" 4 L e 

-l 
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100 km s" 



(2) 



Finally, we can calculate the mass flux through an area A, 
when the gas moves at a velocity v: 



M = pvA — jimnnvA , 



(3) 



where p = fimnn is the mass density, calculated from the mean 
molecular weight fi = 1.24, the mass of the hydrogen atom Mh 
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Table 6. Mass flux in the knots 
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and the atom/ion number density n. The difficulty here is to es- 
timate n from the measured electron density n e . This can be 
done using ionization < I > and the compression < C > from 
the shock models of Hartigan et al. (1994). These authors rec- 
ommend to scale the density with V< C > because this is more 
appropriate for a clumpy medium than scaling with the compres- 
sion factor itself. We get: 



Since we do not measure the size of the jet perpendicular 
to the direction of the flow, we use a radius of 0.15", which 
is the width of the innermost knots resolved by Wassell et al. 
(2006). We take the total space motion of knot B as the flow 
speed of material. Table 6 gives the mass flux in the knots B and 
B2 estimated from all three methods. Given the simplifications 
the difference in the mass flux measured is not surprising. The 
densities of knot B and B2 are very similar, but their luminosities 
differ. This indicates, that we see a larger emitting volume in 
knot B2. In knot B the mass flux calculated from the flow is 
much larger than the mass flux from the line luminosities. This 
could indicate, that either a larger part of knot is too cool to 
radiate in [O i] or [S n] or that the knot has a small filling factor, 
i.e. part of the volume contains cool material unrelated to the jet 
(Podio et al. 2006). 

5. Discussion 

5.1. The central star 

HD 163296 is remarkably constant over nearly 200 ks of X-ray 
observing time. This is consistent with the abundance pattern of 
an inactive star and the cool temperatures. In contrast, CTTS 
typically are active stars with frequent flaring and an inverse 
FIP abundance pattern. A jet collimation shock could explain 
some of the soft emission. However, no offset of the X-ray emis- 
sion from the central star is seen and we estimate an upper limit 
of 30 AU. This is stricter than the number given in Gunther & 
Schmitt (2009) and we can thus calculate a more stringent limit 
on the conditions at the jet base, assuming that a shock at the 
jet base is the origin of the central soft X-ray component. The 
post-shock cooling length d coo i depends only on the pre-shock 
density « an d shock temperature and can be written as (Raga 
et al. 2002): 

j OATT/ 105cnr3 \/ ^hock "I 45 

c/ c00 i*8Au T^n r • (5) 

\ «o /UOOkms- 1 / 



Gunther & Schmitt (2009) derive u sh0 ck = 400 km s _1 , thus the 
lower limit on the density is «o > 2 x 10 4 crrT 3 . The flow time 
through this post-shock region would be about a year. In turn, 
the density sets an upper limit on the radius R of the shock front 
assuming a circular shock at the base of a cylinder (for details 
see Gunther et al. 2009) 

, ZlOWn ' 5 / VEM f 5 ( 0.2 lkeV\ 1125 

where R is given in AU. With the VEM from table 2 this leads to 
an upper limit R < 5 AU. This radius would require a rather large 
opening angle at the base of the jet, thus it is likely that the actual 
densities are higher than the lower limit derived here. The mass 
flux required to power the X-ray emission is 6 x 10-"M o yr -1 , 
only about a tenth of the mass flux in the outer optical knots (ta- 
ble 6). Thus, this seems to be powered by a different (faster and 
less massive) component of the outflow than the optical knots in 
the jet. 

5.2. The one sided Lya emission 

We find a significant extension of the Lya line in the direction of 
the approaching jet in all observations. The signal on the blue- 
shifted side of the emission line could partly be caused by the 
jet itself, but the red-shifted emission has to come from a dif- 
ferent origin. The emission is seen between 0.1" and 1", which 
corresponds to 12 and 120 AU if the emitting structure is in the 
plane of the sky. Perpendicular to the jet the emission is mostly 
contained in the inner 0.2". The emission reaches as far out 
as 400 km s -1 as seen e.g. at the position of the first knot (see 
Fig. 6, where this feature is labeled "scattered emission?") but 
any signal within 150 km s _1 of the line center is hidden behind 
interstellar absorption and geocoronal emission. 

This feature was noticed already by Devine et al. (2000) who 
suggest it to be caused by infall or a poorly collimated wind. We 
disagree with this explanation because of the large width of the 
observed line. Thermal broadening cannot produce the width 
since hydrogen would be fully ionized at the required tempera- 
tures and pressure broadening only becomes relevant for densi- 
ties of the order 10 13 crrT 3 (Vidal et al. 1973). Thus, the line 
width must be explained by kinematic motions but at distances 
of a few AU where the feature is seen the gravity of the central 
star provides an accretion flow speed of a few km s _1 only, much 
less than the width of the emission line. 

The disk obscures the far side of the system, so that we only 
see the approaching wind. If we were looking pole-on all out- 
flows on the visible side would be blue-shifted, but since the 
inclination is 50-60° some outflows might appear red-shifted. In 
the extreme case of a flat disk and a disk-wind moving in par- 
allel to the disk the flow speed still has to reach 500 km s _1 in 
order to achieve a line-of-sight velocity of 400 km s -1 . Not only 
is the limiting case of a flat disk unrealistic, the required ve- 
locity for this outflow is also much larger than the jet velocity. 
Although we do not understand the jet collimation process in de- 
tail, current models suggest that the collimation happens within 
a few stellar radii (e.g. Lii et al. 2012) and commonly the less- 
collimated layers of the outflow are much slower. One mecha- 
nism for such an outflow is a photo-evaporative disk wind, where 
the upper layer of the disk is heated by the stellar radiation. This 
layer is hot enough to emit in Lya, but predicted velocities for 
CTTS are only a few km s _1 (Ercolano & Owen 2010; Pascucci 
et al. 2011). Even for HAeBe, which are much brighter in the 
UV, these winds will not reach the required velocities. Addi- 
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tionally, all the above models would predict an equally strong 
blue-shifted emission which is not observed. 

These difficulties make it promising to propose a model, 
where the Lya flux seen at some distance to the central source 
does not originate there but is scattered emission from the central 
source, since a spectrum extracted at the stellar position shows 
that the Lya line is wide enough. Presumably, most of the hydro- 
gen emission originates in the accretion funnel and the accretion 
shock (Muzerolle et al. 1998). However, dust scattering is not 
wavelength selective and would thus lead to comparable spatial 
profiles for the FUV continuum and the wings of the Lya line in 
contrast to Fig. 7. (The same argument makes an instrumental ar- 
tifact unlikely.) This can be overcome, if the Lya line is formed 
in the accretion funnels above the plane of the disk. The stellar 
continuum would be absorbed by the inner disk wall, while the 
accretion funnels could reach above the rim and illuminate the 
disk surface out to the required radius. Wisniewski et al. (2008) 
observe dark lanes in scattered optical light a few arcseconds 
from the central source, which correlate with the overall disk lu- 
minosity and are consistent with shadowing by a time-variable 
inner disk wall. An alternative explanation is a dusty disk wind, 
with dust grains entrained in the outflow (Bans & Konigl 2012). 
A flared disk or a dusty disk wind would intercept a large frac- 
tion of the total Lya flux, tn principle, we should see Lya from 
the entire disk surface but it is only observed in the negative spa- 
tial direction (below the star in Fig. 5 and to the left in Fig. 7). 
This is the far side of the disk that we see through the jet, so 
the light needs to be back-scattered to reach us. However, in the 
UV dust is strongly forward scattering (Draine 2003) and thus 
the strongest signal should be observed in the positive spatial di- 
rection. The scattering angles in a dusty wind are slightly more 
favorable than in a flared disk. 

Alternatively, scattering by dust in the jet would explain nat- 
urally why we see the extended emission on the side of the jet 
only and also why it seems not to be extended perpendicular to 
the jet axis. Also, the required scattering angle is the inclination 
of the jet (about 50°), which is in the preferred range of angles. 
A potential problem with this scenario is that there is no direct 
observational evidence for dust in the jet. Theoretically, if the 
jet material is launched close enough to the central star, all dust 
may be destroyed near the jet base, by stellar irradiation and/or 
thermal evaporation, although detailed calculations are needed 
to quantify the dust content of the jet. 

Since the spatial profiles from observations between 1999 
and 2011, 12 years apart, are compatible it is unlikely that a 
feature in the disk, e.g. a density enhancement on the far side of 
the disk, causes the scattering to be much more efficient on one 
side than on the other. 

We conclude that several components of dust scattering, 
wind emission and wind absorption with different speeds and 
launching radii need to be present at the same time to produce 
the observed signature. 

5.3. Time evolution of the knots 

Knot B is the only knot observed in 1999 and in 2011. In sec- 
tion 4.3.1 we already derived its proper motion. In this time 
period the observed surface luminosity decreased significantly. 
This could be due to intrinsic dimming or indicate that the jet 
width increased beyond the width of the slit used in the obser- 
vations. The lowest difference -a factor 2- is observed for [Oi], 
the largest difference -a factor 10- is observed for Ha. Conse- 
quently, the ionization fraction and shock speed deduced from 
the line ratios are lower in 201 1 than in 1999. 



Except for the [S n] doublet no element is observed with 
more than one line, thus the temperature of the emitting material 
cannot be measured independent of the abundances. However, 
since all observed elements are neutral or singly ionized, plasma 
with temperatures between 5000 and 25000 K must be present. 
In this energy range, atomic lines are the dominant contributors 
to line cooling. As an estimate, we use the radiative loss rate 
A = 5 x 10~ 23 erg cm 3 s -1 which is valid for this temperature 
range within a factor of a few. The cooling time r can be calcu- 
lated as the ratio of the thermal energy content of the gas and the 
radiative loss as 



T = 



3/2nkr 
Knn e 



0.4 yr , 



(7) 



where k is the Boltzmann constant and we use the n e and n es- 
timated in Sect 4.3.3. The short time scale shows that the knot 
must be continuously reheated by a shock that converts the ki- 
netic energy of the jet into thermal energy. Ultimately, the ra- 
diation is thus powered from the kinetic energy of the jet. The 
kinetic energy density of the material is 300 times larger than 
the thermal energy density. Indeed, v s k y and the shock velocity 
are both lower in the 2011. This is consistent with a jet that de- 
celerates as its kinetic energy is radiated away. An additional 
consequence is that the measured proper motion, which is deter- 
mined from the position of the shock front, is not the same as the 
proper motion of the mass in the jet. 

The knots A and C are not visible in the new observation. 
As shown above the cooling times for gas in the knots are very 
short, so no detection is expected, if the shock fronts in these 
knots are no longer strong enough to power the radiation. 

Due to the low signal in the X-ray observations we cannot 
derive meaningful limits on the density and cooling time, but 
we note that in DG Tau the outer X-ray knots are less luminous 
than the inner ones (Giidel et al. 2008) and a cooling of the X- 
ray plasma as it travels outward has been observed by Schneider 
etal. (2011). 

5.4. Long term evolution of the jet 

In ground-based observations Wassell et al. (2006) are able to 
identify fainter knots further away from the driving source and 
they estimate a dynamical age of HH 409 of at least 80 yr. Our 
observations are not as sensitive to fainter structures, but in the 
inner jet we can confirm that the general shape of the jet has 
not changed over one decade. We have discovered a new knot 
(B2), which has a similar velocity as the previously innermost 
knot B. The distance to the preceding knot is also comparable, 
indicating a periodicity for knot launching. However, knot B2 
has a lower shock velocity and compression factor, but slightly 
larger electron density and mass loss rate than knot B had at a 
comparable distance from the driving source. This shows that 
the launching conditions must have changed slightly. 

5.5. Implications for jet launching 

Most theoretical models rely on magnetic fields to explain the 
jet launching and collimation. Yet, it is unclear if the central star 
HD 163296 posses a large scale magnetic field. On the one hand, 
none was found in the observations of Hubrig et al. (2007), on 
the other hand, we observe a constant X-ray flux with tempera- 
tures above 2 keV, which is beyond the heating from an accretion 
shock and is presumably caused by magnetic reconnection and 
thus indicates the consistent presence of stellar magnetic fields. 
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We know that knot B2 was launched within a few years of the 
optical polarimitry and the X-ray observations: If it has moved 
with constant proper motion, then it was launched in 2002, while 
the observations of Hubrig et al. (2007) were performed in 2004 
and we analyzed X-ray data from 2003, 2007 and 2011. If the 
limit on the magnetic field is representative then a strong stellar 
magnetic field is not a necessary pre-requisite of jet launching. 
Instead, Hubrig et al. (2007) report the detection of Stokes com- 
ponents in the Ca H and K lines, which are much wider than ex- 
pected from the photosphere. If these lines are of circumstellar 
origin they probe the jet launching and collimation region more 
directly. These fields could be generated in e.g. a disk dynamo 
(Vishniac & Brandenburg 1997). 

Alternatively, we can take the presence of hard X-rays as 
proof of stellar magnetic fields. The optical polarimitry concen- 
trated on hydrogen lines, which - at least in CTTS - have mul- 
tiple contributions from the inner edge of the disk, the accretion 
funnel and the star. The stellar contribution of the flux might be 
too small to produce a measurable Stokes signal. Thus, we argue 
that the jet launching and collimation in HAeBes is comparable 
with the process in CTTS. 

6. Summary 

We analyze a series of Chandra and HST observations of 
HH 409, the jet from the Herbig Ae star HD 163296. In X-rays 
the star itself is very stable and the observed temperatures and 
emission measures are fully compatible with previous observa- 
tions. No extension beyond the PSF is found for the central star, 
which limits any unresolved sources with a significant fraction 
of the total luminosity to a maximum separation of 30 AU. This 
is important, because an XMM-Newton spectrum showed that the 
soft X-ray emission originates above the stellar surface, possibly 
in the jet collimation region. While in 2003 a weak feature on 
the edge of the stellar PSF is observed at a distance of -4", it is 
unclear if emission at that position is present in 2011. Swartz 
et al. (2005) claim the detection of a knot in the jet at -7.2"; 
our analysis confirms this. There is a tentative grouping of three 
photons at -10.6" in 2011. This position is consistent with the 
2003 knot, if it propagates with the same proper motion as the 
optical knots. In any case, the luminosity of the knot dropped 
significantly since 2003, indicating a short cooling time. 

In the UV the Lya line is the most prominent feature. In 
the new observations we detect two knots in the jet, but not the 
counter-jet. In addition, there is extended emission on the side of 
the approaching jet that is red-shifted up to 400 km s _1 and can 
be seen out to about 1". On the blue-shifted side the emission 
is much weaker. Neither red nor blue-shifted emission is seen in 
the direction of the counter-jet. It requires multiple components 
of accretion, a wide-angle wind or dust scattering of the Lye line 
from the accretion spot to explain this feature. 

In the optical we detect Ha and several forbidden emission 
lines. We detect a new knot, called B2, and re-observe the previ- 
ously detected knot B. The distance between those two knots is 
similar to the distance between knot B and the preceding knot. 
This nicely fits into a series of knots in the counter-jet with al- 
most equal spacing which were launched over several decades 
(Wassell et al. 2006). Knot B moves with a proper motion com- 
patible to the previous values for knot C. Knot B is significantly 
fainter now than in 1999, particularly in Ha. The line ratios indi- 
cate a lower temperature. A significant spatial expansion, which 
would bring the size of the knot beyond the width of the slit, 
could also contribute to the fading. Still, the cooling time scale 
is less than a year, so knot B needs to be continuously re-heated. 



This supports a picture where the knots are shock fronts traveling 
though a jet of continuously launched, cool material. Over time, 
these shocks radiate away their energy, so the post-shock cooling 
zone formed behind them has consecutively lower temperatures. 

Knot B2 today is seen at almost the same distance from the 
central star as knot B was a decade ago. Yet, knot B2 has a lower 
ionization and radial velocity, but a larger mass flux than knot B 
had a decade ago. Thus, the jet launching conditions must have 
changed in this period. 
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wavelength [A] 

Fig. 1. G750M data from 201 1. Hot pixels close lines of interest have been removed manually (see text for details). Lines which we find to be 
spatially extended (see Sect. 4.3) are marked. 




wavelength [A] 

Fig. 2. G140M data from 201 1. The Lya line airglow feature has been removed. Lines which where found to be extended in earlier observations 
are marked. 



